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Abstract. In a recent paper we introduced a multivariable generalisation of 
the Bessel polynomials, depending on one extra parameter, and related to the 
so-called hyperbolic Sutherland model with external Morse potential. In this 
paper we obtain a corresponding multivariable generalisation of a well-known 
orthogonality relation for the (one-variable) Bessel polynomials due to Krall 
and Frink. 



1. Introduction 

Up to normalisation, the Bessel polynomials ym{x', a), which depend on a (com- 
plex) parameter a and a non-negative integer m, can be defined as the polynomials 
of degree m such that they are eigenfunctions of the ordinary differential operator 

As such, they first appeared in a paper by Bochner [Boc29], although the name 
Bessel polynomials was introduced later by Krall and Frink [KF49] . We recall that a 
second parameter b can be introduced by a change of coordinates to y = bx. As first 
shown by Hahn [Hah35], there exist no interval on the real line and corresponding 
non-negative weight function with respect to which the Bessel polynomials are 
orthogonal. However, Krall [Kra41] demonstrated that the Bessel polynomials are 
orthogonal in a more general sense. In fact, they satisfy the orthogonality relation 

(1) / yn{x)ym{x)p{x)dx = 0, n ^ m, 

J\x\ = l 

where 

r(a) / 2 



27rj ^-^ T(a + n — 1) V x 

n—O ^ ' ^ 

The result was first stated in this form by Krall and Frink [KF49]. For a general 
account of the theory of Bessel polynomials, and an excellent guide to the extensive 
literature on the topic (up to 1978), see the book by Grosswald [Gro78]. 
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In a recent paper we introduced and studied a multivariable generalisation of 
the Bessel polynomials as eigenfunctions of the partial differential operator 



dijC: dx^ Xf^ X 't dxf^ 

where k is a (complex) parameter; see Hallnas [Hal08]. In this latter paper we 
emphasised the relation with a quantum many-body system of so-called Calogero- 
Moser-Sutherland type. This lead to certain orthogonality results with respect to 
the inner product in a particular Hilbert space of square-integrable functions, to 
be discussed further in Section 13.31 However, these orthogonality results are valid 
only for a non-generic set of parameter values, and then involve only a finite set 
of multivariable Bessel polynomials. The main purpose of this paper is to deduce 
a multivariable generalisation of the orthogonality relation ([T]), valid for generic 
parameter values, and involving all multivariable Bessel polynomials. 

We conclude this introduction by mentioning that Duran [Dur93], Evans et 
al. [EEKL93], as well as Kwon et al. [KKH92], have obtained real, but necessarily 
not non-negative, weight functions with respect to which the (one-variable) Bessel 
polynomials are orthogonal on [0,c»). The problem of obtaining corresponding 
multivariable weight functions remains open. 

2. Orthogonality relation and normalisation factors 

In this section we state our main result: the orthogonality relation for our mul- 
tivariable Bessel polynomials. In order to do so we first briefly review some basic 
results on symmetric functions and Jack polynomials, recall the definition of our 
multivariable Bessel polynomials, and state a result on the nature of their depen- 
dence on the two parameters a and k. 

2.1. Symmetric functions and Jack polynomials. A partition A = (Ai, . . . , A„) 
is a sequence of non-negative integers such that Ai > • • • > A„, and its weight is 
the sum |A| = Ai -I- • • • -I- A„ of its parts. On the set of partitions of a given fixed 
weight the so-called dominance order is defined by 

< A <^ ^1 + • • • + /ij < Al + • • • + Ai, Vi = 1, . . . , n. 

The monomial symmetric polynomials m\ in the variables x — (xi, . . . ,Xn) are 
given by 

■mx{x) ^^x"' ■ • 

a 

where the sum is over all distinct permutations a of A. These monomials form a 
linear basis for the algebra of symmetric polynomials in n variables with complex 
coefficients, from hereon denoted A„. As long as the parameter k is not a negative 
rational number, there exists for each partition A — (Ai, . . . , A„) a unique symmetric 
polynomial P\ in the variables x, referred to as a Jack polynomial, of the form 

P\ = mx + ^ uxf^nifj, 

for some coefficients ux^ such that it is an cigcnfunction of the partial differential 
operator 



D=y xi—^^2K 



dx- . . Xi Xj dxi 
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see e.g. Section 3 in Stanley [Sta89]. We remark that the parameter k is the inverse 
of the parameter a used by Stanley. We let T" = {x G C" : |xi| = l,Vi = 
1, . . . ,n} be the n-dimensional torus. In addition, we will use the notation = 
[l/xi, . . . , 1/xn), dx = dx\ ■ ■ ■ dxn and dx/x = dx\/xi ■ ■ ■ dxn/xn- If > then 
the Jack polynomials are orthogonal polynomials with respect to the bilinear form 
on An defined by 

for any f^g^ An-, and the corresponding normalisation factors are given by 
(P py ^ n r(Ac(j - z + 1) + A, - Xj)TiKij - i - 1) + A, - + 1) 
i<y<n - + Ai - A,)r(K(i -i) + Xi- Xj + 1) ' 

see e.g. Section VI. 10 in Macdonald [Mac95]. Here, as well as in the discussions 

below, we use the convention that empty products arc equal to one. For a detailed 
exposition of the theory of symmetric functions and Jack polynomials we refer to 
the book by Macdonald [Mac95]. 

2.2. The mult i variable Bessel polynomials. In analogy with the Jack polyno- 
mials, the multivariable Bessel polynomials can now be defined by the following 
fact (see Section 3 in Hallnas [Hal08]): given a partition A = (Ai,...,A„) there 
exists for generic values of the parameters a and k a unique symmetric polynomial 
Yx of the form 

for some coefficients ux^^ such that it is an eigenfunction of . We remark that it 
is rather easy to establish the existence of such a polynomial for generic parameter 
values, but to obtain a condition on the parameters which is both necessary and 
sufficient is much harder. In Hallnas (sec Proposition 3.1) we obtained a sufficient 
condition which, however, was rather far from being necessary. To state our main 
result we will need such a condition which is much closer to being necessary. We 
will resolve this problem by proving the following: 

Proposition 2.1. The multivariable Bessel polynomials Yx{a,K) are rational func- 
tions of a and k with poles only along hyperplanes given by k a negative rational 
number or a — 1 + ki a non-positive integer for som,e ?' = 0, 1, . . . , 2(ri — 1). 

2.3. Statement of main result. We will show that the multivariable Bessel poly- 
nomials are orthogonal polynomials on the torus T" with respect to the weight 
function 

Xi 



W{x;a,n) = 'W{x;a,n)'[\_il 



where 

m T(x-a;z) ^ Vfn n<3-i) + Xi-X, + l) 



{2mYn\ ^ \ 11 r(K(i - i - 1) + A, - A,- + 1) 



A r(^^)r(a + ^c(2n-j-l)) \ 1. 
y r(K)r(a - 1 + K(2n - i - 1) + Ai) J ' ^' 
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A priori, it is not clear for which x e C" the right-hand side of ([2]) converges, or 
indeed if such x exist. This problem is settled by the following: 

Lemma 2.1. Assume that k is not a negative rational number and that a — 1 + 
K(2n — i — 1) is not a negative integer for any i = 1, . . . ,n. Let N be a neighbourhood 
of € C" . Then the right-hand side of ^ converges uniformly and absolutely on 
C" \iV. 

We let [a\m denote the Pochhammer symbol, defined by [a\o = 1 and [a\m — 
a{a + 1) • • • (a + m — 1) for all positive integers m. In addition, we will use the 
notation [ai, . . . , ap]„i — [ai]m ■ ■ ■ [ap]m- We are now ready to state our main 
result. 

Theorem 2.1. Assume that 5R(k) > and that a — 1 + ni Q, —1, . . . for all 
i = 0, . . . , 2(n — 1). Let A = (Ai, . . . , A„) and \x — (/ii, . . . , be two partitions 
such that X ^ fj,. Then 



Remark 2.1. Suppose that n ^ I. Then, as is readily inferred from their definitions, 
P{m)ixi) = and Y(m)(a;i) = Cmym{xi) for some constant Cm, depending on the 
particular normalisation chosen for the ym- We recall that Krall and Frink used 
the normalisation in which the constant term of each ym is equal to one. In this 
case. Cm = 2'"[a — llm/l^ — l]2m; see Section 6 in Hallnas [Hal08] for the precise 
relation between this and the monic normalisation in the multivariable case. It is 
thus clear that W{xi) — p{xi), and that we recover Krall and Frink's orthogonality 
relation ([T]) for the Bessel polynomials ym by setting n = 1 in Theorem l2.1l 



In this section we prove the results stated in Section [21 

3.1. Proof of Proposition 12.11 In a first step we will deduce a closed expres- 
sion for the coefficients in the expansion of the multivariable Bessel polynomials 
in Jack polynomials. Our starting point will be such an expression for the Jacobi 
polynomials associated with the root system BC„, obtained by Okounkov and 01- 
shanski [OO06] (see Proposition 2.3). The reason is that the multivariable Bessel 
polynomials can be obtained from these BCn Jacobi polynomials via a limit transi- 
tion; see Proposition 4.1 in Hallnas. By applying this limit transition to Okounkov 



(3) 




Moreover, 




3. Proofs 
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and Olshanski's result we will obtain the desired expression for the multivariable 
Bessel polynomials. It will then be a straightforward exercise to infer Proposition 
12.11 from this expression. 

We let t{z) = (4sinh^ zi/2, . . . ,4sinh^ •z„/2). For each partition A — (Ai, . . . , A„) 
and generic (complex) values of the parameters (fci , fe, fca) a BCn Jacobi polynomial 
P^'^' can be defined as the unique eigenfunction of the partial differential operator 

+ fc3 g (coth - z,) (|- - ^) + coth + z,) (|- + ^)) 
such that 

P^''{z) = P,(t{z)) + Y,u,Mt{z)) 

for some coefBcients w^^; see e.g. Beerends and Opdam [B093] and references 
therein. An important ingredient in Okounkov and Olshanski's result is the so- 
called BCn interpolation polynomials. We proceed to recall a 'combinatorial' for- 
mula for these polynomials which we shall make use of at a later stage in the proof. 
To this end we let fi = (fii, . . . , be a partition, and recall that a reverse column- 
strict tableux T of shape fi can be defined as a function that assigns to each box 
ihj) G an integer T(i,j) G {1, . . . ,n} such that T[i,j) decreases strictly down 
each column and weakly from left to right along each row. Suppose that k is not 
a negative rational number. The BCn interpolation polynomial can then be 
defined by a formula of the form 

(5) /i) = ^?At(k) n {^^T(t,j) +h- HiT{i,j)f 

~ {j ^ I - K{i - \) + h - KT{i,3)f), 

where the sum is over all reverse column-strict tableux T of shape /j,; see Proposition 
2.2 in Okounkov and Olshanski [OO06]. The weights iPt{k) are rational functions 
in K with poles located only at negative rational numbers. An explicit formula for 
these weights can be found in Section VI. 10 of Macdonald [Mac95] (see Equations 
(10.10)-(10.12)). Since we will not make use of this formula we refrain from stating 
it here. The results of Proposition 2.3 in Okounkov and Olshanski [OO06] can now 
be formulated as follows: 

Theorem 3.1 (Okounkov and Olshanski). Let a = {ki + 2fc2)/2. Then 

Pf^(z; fci,fc2,fc3) ^ /^(A; fca, cr + k3n)Pf,{t{z); fca) 

^ ' Pf^(0";fci,fe,A:3) ^^I^{^l;k3,a + k3n)pBC(on.k^^k^^k3)■ 

Remark 3.1. Okounkov and Olshanski referred to this formula for the BCn Jacobi 
polynomials as a binomial formula. 

At this point we set (fci, /c2, ^3) = ((a-l-|-2e^)/2, (a-l-2e')/4, k) with e > 0. In 
particular, we then have cr = (a — l)/2. We have, furthermore, the limit transitions 

lim e-"l^lpf ^(zi +e,...,zn + e; fci, /c2, ^13) = Yx(e'' , . . . , e^" ; a, k) 

e— >oo 



6 



MARTIN HALLNAS 



and 

lim e-^l^lpf ^(0"; h,k2, h) = yA(0"; a, k); 

e — >QQ 

see Proposition 4.1 and the proof of Proposition 6.1 in Hallnas [Hal08], respectively. 
We observe that hnig^^oo 46"*^ sinh'^(z + e)/2 = e^. By taking the hmit e ^ cx) in 
^ we thus obtain the foUowing: 

Proposition 3.1. The multivariable Bessel polynomials have the series expansion 
Yx{x;a,K) y-^ I ^{X; k, (a ~ l)/2 + Kn)Pf^{x] k) 



_ . J- II, 



Yx{0^; a, k) ^ /^(^; (a - l)/2 + K,n)Y^{0"; a, k) ' 

We recall that the Jack polynomials Px{k) are rational functions of k with poles 
only at negative rational numbers; see e.g. Equation (10.13) and the preceding 
discussion in Chapter VI of Macdonald [Mac95] . We observe that (O implies that 
this fact holds true also for /^(A; k, (a — l)/2 + ku). As discussed by Okounkov and 
Olshanski [OO06], the BC„ interpolation polynomials are normalised such that 

I^ip; K, {a - l)/2 + nn) 

= (l + /ii - j + K{fi'j - i)){a - 2 + fii + j + K{2n - ^'^ - i)) . 

It is clear that a factor 1 + /ij — j + nip-'j — i) = only if k is a negative rational 
number. Since /i^ + j — 1 > and 2 < ji'^ + i < 2n, it is also clear that a factor 
a — 2 + fii+j + K{2n — fj,j — i) = only if a — 1 + is a negative integer for some 
i = 0,...,2{n — 1). Furthermore, the specialisation of the multivariable Bessel 
polynomials at a; = (0") is given by 

Y m"- n k) - n - ^ + 1)]a.-a, K2n -1-3 + 1)]a.+a,- 



n 



[a — 1 + K(n - i)]\^ 



[a - 1 + 2K(n - i)]2\i ' 

see Proposition 6.1 in Hallnas. It is thus clear that the zeros, as well as the poles, 
of Ya(0"; a, k) are located along hyperplanes given by k a negative rational number 
or a — I + Ki equal to a non-positive integer for some i = 0,...,2(n — 1). This 
concludes the proof of Proposition l2.1l 

3.2. Proof of Lemma 12.11 In order to prove the statement we will make use of 
the so-called integral form 

Jx{x; k) = K'^^^h^Pxix; k), = Y[ Y[{K - 3 + i^{Xj - i + 1)), 

j=i j=i 

of the Jack polynomials. We recall that the coefficients in their expansion in 
the monomials m\ are polynomials in 1/k with positive integer coefficients; see 
e.g. Knop and Sahi [KS97]. Since the Jack polynomials J\ are homogeneous of 
degree |A|, it follows that 

|Ja(x;«)1 < Ja(1";|k|)|N|I^I 
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with \ \x\\ = maxQxiDf^i. We recall from Stanley [Sta89] (see Theorem 5.4) that 

A. 

Ja(1"; n) = - I + tiin - i + 1)). 

We observe that since k is not a negative rational number, there exist e > such 
that \Ki + m\ > e for all i = 1, . . . , n — 1 and positive integers m. We fix such a 
number e. ft is then clear that 



n 

i<j 



Kj + 1]a.-a, 



-n n 

i<ij rrtij—l 



- i) + "f^i] 

Ai-Ai 



< R 



A| 



for some i?i independent of A. Proceeding in a similar fashion, and using the fact 
that J(l"-; \k\) < [1^=1 II'^I^Ia,, it is readily verified that also 

l^l'^'JA(l";|Ai|) , p|A| 

m - ' 

for some i?2 independent of A. We thus conclude that for some constants C and R, 
independent of A, we have 



\W{x)\<CY, 



r(a + K(2n-2- 1)) 



A 

n oo 



< 



1=1 Ai=0 



n r(a - 1 + K(2n -i-l) + X,) 



|r(a + K(2n- j- 1))| 
|r(a - 1 + K(2n - i - 1) + A,) 



which clearly converges uniformly and absolutely on C \ iV. The statement thus 
follows from the so-called Weierstrass M-test; see e.g. Theorem 11.5 in Apostol 
[Apo67]. 

3.3. Proof of Theorem 12. li Let v — (y\^ . . . ,Vn) be a partition. Using the 
orthogonality of the Jack polynomials, as well as the fact that x^ ■ ■ ■ x^Pv^^x) — 
P,y+(in) and Py is homogeneous of degree \v\^ it is a straightforward exercise to 
verify that 



(7) 



Py[x)W{x)dx 



r(K(j - i + 1) + Ml - Vj) 



I2\\v\+n TT , 

^ ' V{K[j -i) + v,-Vj) J^-J- r(K)r(a + K(2n - i - 1) + i/j) ■ 



n 



T{Ki)T{a + K{2n-i - 1)) 



Since T{a + m)/T{m) ~ [a]™, it follows from ([7]) and Proposition 12. II that 



(8) 



Y^{x)Y^{x)W{x)dx 



n 



T{Ki)T{K{i - i)) 



\f{a.,K) 



for some rational function / with poles only at the excluded values of a and k. 
In order to verify that the right-hand side of ([H]) is zero, i.e., to establish the 
orthogonality relation ([3]) , it is thus sufhcient to verify that this function / vanish 
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on a subset of which is dense in the Zariski topologyQ To this end we fix a non- 
negative integer m such that m > max(|A|, and consider the set of parameter 
values (a, k) G such that k > and a < — 2(m + K{n — 1)) + f. It is clear 
that this subset of is dense in the Zariski topology. We recall from Theorem 7.1 
in Hallnas [Hal08] the following orthogonality relation for the multivariable Bessel 
polynomials of degree at most m: if A and are two partitions such that |A|, < m 
and \ ^ iJi then 



(9) / Y^{x)Y^{x)WL2{x)dx = 0, WL2{x) = -\{xr^e-^'-^\{\x^,~x,\^\ 

1=1 t<j 

We stress that if either of the partitions A and /i has a weight greater than m then 
this integral might not exist; see Lemma 7.1 in Hallnas. We proceed to compute the 
integral Pu{x)Wi,2{x)dx for the partitions v — {vi, . . . , i/„) such that < 2m, 
and to show that, up to a non-zero constant, it is also given by the right-hand side 
of ([7]). As we will show below, the value of this integral is easily inferred from a 
particular limiting case of the following generalisation of Selberg's integral formula 
due to KadeU [Kad97]: 

Theorem 3.2 (KadeU). Let v = (z/i, . . . , he a partition and let a, /?, k £ C be 

such that 



(a) > -Vn, > 0, > 0. 



Then 



1 r " 

(10) - / PAy;^)l[y^~'(^~y^^^~']l 

"'[0.1]" i=l i<j 



_ T-r r(K(j - i + I) + i^i ~ Vj) -A- T{a + K{n - i) + + K{i - 1)) 

- 11 r(K(j + i^j) l\ T{a + P + K{2n-i-l) + u,) ' 

We recall that (1 — x/ 0)^ — > as (3 ^ oo, and also that 

H-^oo' r(/3 + 6) 

for any complex a and 6; see e.g. Theorem 3.4-1 in Carlson [Car77] for the latter 
statement. By first replacing each variable yt by ?/j//3, then using the fact that 
Pu is homogenous of degree |^|, and finally taking the limit [3 — s- oo, it is thus 
straightforward to infer from PU)) that 



1 r " 



n T{K{j - i + 1) + - t^j) -/V , , ^ 

r(«(,-.) + ..-^) n r(« + <n - .) + .0; 



^We recall that a subset y C is dense in tlie Zariski topology if and only if the zero 
polynomial is the only polynomial which vanish on Y. 
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c.f. Section 17.3 in Mehta [Meh91]. On the other hand, by a change of variables to 
(2/1, . . . , ?/„) = (2/xi, . . . , 2/a;„), it is easy to verify that 



n{a—l)n-{-K,n{n—l)-\-\i/\ r 



i<2 

We fix a non-negative integer TV such that v C (A^"), and let ii be the partition 
obtained by taking the complement of v in (iV"). Explicitly, 

i>i - i = l,...,n. 

Then 

Indeed, the two sides have the same leading term, and it is readily verified that they 
satisfy the same orthogonality relation; see Macdonald [Mac]. By combining the 
observations above it is now a straightforward exercise to establish the following: 

Proposition 3.2. Let v ~ (i^i, . . . , Vn) be a partition and a, k £ <C be such that 

3?(k)>0, K(a) < -i/i -23fi(K)(n- 1) + 1. 

Then 

I P.{x)Wl. {X)dx = 2^-~^)n+.n(n-l) + H TT T - t + I) + - l^j) 

Jr^ r{K{j - 1) + Ui - i^j) 

n 

Yl r(-a + 1 - K{2n - i - 1) - t/,). 



X 

1=1 



We observe the formula T{~a - m + 1) = (-1)"T(-q; + l)r(a)/r(Q; + m). By 
comparing Proposition 13. 21 with ([7]) we thus find that 



(11) / P^{x)W{x)dx 



T" 



2ia-2)n+Kn{n-l) 11 Y{Ki)r{-a + 1 - k(271 - i - 1)) jRr^ 



if the conditions on a and k stated in Proposition 13.21 are satisfied. We observe 
that the product of gamma functions preceding the integral in the right-hand side 
of pTjl is non-zero. Since the Jack polynomials P^, with < 2m span the space of 
symmetric polynomials of degree at most 2m, the orthogonality relation ^ thus 
implies that the function / in the right-hand side of ([5]) vanish for all (a, k) G 
such that K > and a < — 2(m + K{n — 1)) + 1, and hence is identically zero. 
This concludes the proof of the orthogonality relation ([3]). Moreover, by using (fTTjl . 
and proceeding as above, it is now straightforward to infer the normalisation factor 
formula (j4]) from Lemmas 7.3 and 7.4 in Hallnas [Hal08]. 
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